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Introduction 

The  present  phase  of  the  research  project  concentrated  on  quantitative 
evaluation  of  the  ceramic  particles  induced  into  the  porous  substrate  as  a 
function  of  deposition  and  solvent  parameters.  This  evaluation  was  carried  out 
both  experimentally  and  theoretically  based  on  an  analytical  model  for  deep 
electrophoretic  penetration  of  particles  into  porous  substrates.  In  parallel 
to  the  electrophoretic  studies,  the  electrolytic  deposition  of  ceramic 
coatings  from  aqueous  solutions  was  studied  from  the  electrochemical  and 
microstructural  aspects.  The  four  enclosed  manuscripts  describe  in  detail  the 
above  mentioned  studies.  Following,  therefore,  is  a  concise  description  of 
the  studies  and  their  results. 

Experimental 

The  experimental  studies  were  carried  out  on  porous  graphite  substrates  - 
UCAR  Grade  45.  This  substrate  has  a  48%  volume  porosity  with  an  average  pore 
size  of  60/j.  The  ceramic  particles  deposited  and  induced  by  electrophoresis 
were  submicron  colloidal  Si02  (Pyrogenic  Aerosil)  with  a  specific  area  of  450 
m2/g.  The  concentrations  of  the  ceramic  material  in  the  solvent  were  3,  13 
and  30  g //.  The  solvents  in  most  experiments  were  water  and  isopropanol.  For 
the  study  of  the  effect  of  solvent  properties  on  penetration  ethanol  and 
pentanol  were  also  used.  The  field  intensity  was  varied  between  0-300v/cm  and 
the  deposition  time  between  5-120  min. 

For  quantitative  determination  of  the  impregnated  Si02,  the  weight  change 
of  the  specimen  was  recorded,  after  removal  of  the  external  deposit.  The 
amount  of  induced  SiC>2  was  also  determined,  in  some  cases,  by  weighing  the 
residue  after  burn-off  of  the  graphite  at  900°C  for  20  hrs. 
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The  electrolytic  deposition  studies  concentrated  on  the  deposition  of 
Zrt>2  from  a  0.1M  ZrO(N03)2nH2O  aqueous  solution.  The  deposits  were  obtained 
at  current  densities  ranging  from  15-100  mA/cm2  and  durations  of  10-60  min. 

For  kinetic  studies  the  coating  weights  were  determined  with  an  accuracy  of 
0.05  mg.  During  deposition  the  variation  of  cell  voltage,  pH  and  temperatur- 
of  the  solution  were  measured. 

The  microstructure  of  the  Zr02  coatings  was  characterized  by  optical  and 
scanning  electron  microscopy  and  the  phase  content  was  determined  by  X-ray 
diffraction.  The  microstructural  studies  were  carried  out  on  deposits  that 
were  fired  at  400,  600  and  900°C  for  1  hr. 

Theoretical 

The  theoretical  analysis,  at  this  stage,  relates  to  a  single  particle 
travelling  in  a  pore  described  as  a  long  circular  tube  of  mean  pore  diameter. 
Three  driving  mechanisms  are  considered:  the  hydrodynamic  drag  force  exerted 
on  the  particle  due  to  the  electroosmotic  flow  of  the  solvent  inside  the  pore, 
the  electrophoretic  force  exerted  on  the  particle  and  the  stochastic  Brownian 
force  due  to  thermal  fluctuations  of  the  solvent  molecules. 

Summary  of  Results 

In  addition  to  the  formation  of  surface  deposits,  electrophoretic 
induction  of  ceramic  particles  into  the  porous  substrate  was  demonstrated. 

It  was  found  that  solvent  properties  such  as  the  ratio  between  the 
dielectric  constant  and  the  viscosity  (e/?)  have  a  significant  effect  on 
penetration.  Thus  the  amount  of  induced  SiC>2  is  ten  times  larger  for  water 
(e/7  =  81)  than  for  pentanol  ( t/t\  -  420).  Penetration  increases  with  particle 
concentration  and  with  deposition  time.  However,  a  plateau  is  reached  for  the 
dependence  of  induced  Si02  on  deposition  time  due 'to  the  build-up  of  an 
external  coating.  Thus  the  extent  of  penetration  could  be  enhanced  if  the 
external  deposit  were  constantly  removed  or  prevented.  An  optimal  field 
strength  exists  for  particle  penetration  due  to  its  dual  effect.  Increase  in 
field  strength  increases  the  particle  velocity  but  it  also  enhances  the 
build-up  of  the  external  deposit  which  blocks  penetration. 

The  theoretical  analysis  enables  to  predict  the  penetration  depth  of  a 
ceramic  particle  inside  the  porous  substrate  under  the  influence  of  the 
electric  potential  gradient.  It  also  enabled  to  derive  the  non-dimensional 
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parameters  that  affect  the  particle  notion.  It  was  thus  deduced  that  the 
parameters  which  affect  most  the  particle  penetration  are  the  Feclet  number 
(Pe)  and  the  deposition  rate  (1): 

„  Ub  /U'  +  U")b, 

Pe  =  r  =  < — 5 — 5 


where: 

U'  -  electrophoretic  velocity 
U"  -  electroosmotic  velocity 

b  -  pore  mean  radius 

D  -  diffusion  coefficient 

k  -  local  deposition  rate 


Large  Peclet  numbers  and  small  deposition  rates  enhance  penetration.  Thus, 
deep  penetration  (over  hundred  pore  diameters)  can  be  achieved  for  very  low 
deposition  rates  (1  <  0.1)  and  for  very  high  Peclet  numbers  (Pe  >  100). 

Deepest  penetration  is  obtained  for  1=0  which  can  be  achieved  if  a  repelling 
force  between  particles  and  pore  walls  would  exist.  Electroosmosis  increases 
the  Peclet  number  and  thus  penetration.  The  diffusion  coefficient  of  the 
particle  has  a  dual  effect:  its  decrease  increases  the  Peclet  number,  but  it 
also  increases  the  deposition  rate.  The  theoretical  analysis  also  indicates 
how  the  penetration  can  be  controlled  by  selection  of  solvent  properties, 
particle  size  and  concentration  and  electric  field  strength. 

The  electrolytic  deposition  of  Zr02  is  obtained  through  the  following 
reactions: 

a)  dissociation  of  the  zirconyl  salt: 

ZrO(N03)2  -*  ZrO+  +  +  2  N03" 

b)  hydrolysis  of  the  zirconyl  ion: 

ZrO++  +  H20  -s  Zr(0H)2t+ 

c)  interaction  of  the  hydrated  cation  with  OH"  ions  generated  at  the 
cathode  by  reduction  reactions  described  further  on: 

Zr(0H)2+t  +  20H"  -*  Zr(0H)< 

d)  dehydration  of  the  hydroxide: 

Zr(0H)<  -*  Zr02  +  2H20 
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Several  cathodic  reactions  are  possible;  however,  the  cathodic 
polarization  curve  indicates  that  the  reduction  of  NO3'  is  the  predominant  OH" 
producing  reaction.  The  rate  of  deposit  formation  increases  with  current 
density  the  dependence  reflecting  the  rate  of  OH"  generation.  Faradaic 
efficiencies  of  40-80%  only  were  obtained  attributable  to  reduction  reactions 
that  do  not  produce  OH"  (such  as  reduction  of  H+)  as  well  as  to  the  formation 
of  the  hydroxide  at  sites  removed  from  the  substrate.  Cell  voltage  and 
solution  temperature  increase  as  the  coating  process  progresses.  This  is  due 
to  the  formation  of  the  deposit  characterized  by  a  high  electric  resistance. 

The  initial  coating  is  a  Zirconium  hydroxide  gel  with  an  amorphous 
character.  Microcracks  develop  in  the  coating  during  drying  due  to  tensile 
stresses  associated  with  non-uniform  contraction.  Firing  of  the  coating  led 
to  its  crystallization  into  fine,  equiaxed  submicron  particles  of  monoclinic 
and  tetragonal  polymorphs  of  zirconia.  Calculation  of  crystallite  sizes  by 
reflection  broadening,  indicated  them  to  be  2  nm  after  treatment  at  400°C  and 
12  nm  when  fired  at  600°C  for  1  hr.  Under  these  conditions  the  polymorphs 
were  ^distinguishable.  Firing  for  1  hr  at  900°C  resulted  in  distinct 
crystallites  of  both  phases,  20-25  nm  in  size.  The  volume  fractions  were  15% 
and  85%  for  the  tetragonal  and  monoclinic  phases  respectively. 

Future  electrophoretic  studies  will  include  further  optimization  of 
penetration  by  inhibition  of  surface  deposition.  One  approach  will  be  the 
application  of  a  low  cathodic  bias  voltage  to  the  specimen. 

Following  the  derivation  of  optimal  process  conditions  for  penetration, 
experiments  with  lower  expansion  ceramics  such  as  HfTiO^  will  be  performed. 

Further  studies  on  the  electrolytic  deposition  will  concentrate  on 
prevention  of  cracking  of  the  surface  deposit  during  drying  and  study  and 
optimization  of  in-pore  deposition. 
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Abstract 


The  objective  of  the  present  study  was  to  evaluate, 
experimentally,  the  feasibility  of  depositing  ceramic  particles 
on  the  surface  of  a  porous  substrate  and  inducing  them  into  it 
by  electrophoresis. 

It  was  demonstrated  that  it  is  possible  to  charge  and 
electrophoretically  deposit  both  oxide  (fused  and  colloidal 
Si02)  and  non-oxide  (SiC  and  SiN)  ceramics  on  a  conductive 
graphite  substrate.  It  was  also  shown  that  ceramic  particles  are 
induced  into  porous  graphite  and  that  under  optimal  conditions 
the  whole  crossection  of  the  specimens  was  penetrated. 

The  amount  of  induced  material  was  found  to  increase  with  the 
ratio  of  dielectric  constant  to  viscosity(e/n)  of  the  solvent  as 
well  as  with  particle  concentration  and  field  intensity.  However, 
due  to  simultaneous  build-up  of  an  overlying  deposit,  penetration 
as  function  of  deposition  time  reaches  a  plateau  prior  to  filling 
of  the  whole  pore  volume. 
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Introduction 

Electrophoretic  deposition  is  obtained  via  migration  of 
non-conductive  but  electrically  charged  particles,  in  an  electric 
field  towards  one  of  the  electrodes. 

Tvo  different  charging  modes  are  attributed  to  ceramic 
particlesf 1 ] : 

a)  a  dissociation  mode,  such  as  dissociation  of  silanol  groups 
present  in  silica: 

SiOH  +  OH- - >  =SiO-  ~  h20  (negatively  charged  particles) 

b)  an  absorption  mode,  such  as  absorption  of  H*  from  water  or 
from  weak  acids: 

nAl  (OH)  3+HN03 - >(Al(OH)3]n"xH--fxNO-3 

(positively  charged  particles). 

Usually  the  polarity  of  the  particles  has  to  be  determined 
experimentally  as  the  charge  is  acquired  spontaneously  on  mixing 
the  particles  with  the  solvent  and  it  may  be  reversed  upon  addition 
of  ionic  compounds. 

The  surface  charge,  once  formed,  is  then  balanced  by  ions  of  an 
opposite  charqe  derived  from  the  solution  and  thus  a  double  layer 
is  formed  around  it. 

According  to  the  model  proposed  by  Stern,  this  double  layer 
consists  of  a  rigid  part  with  a  linear  potential  gradient,  and  a 
diffuse  part  with  a  non-linear  potential  gradient  termed  "Zeta 
potential "  ( 2 ) . 

When  an  electric  field  is  applied  to  a  charged  particle  the  latter 
tends  to  migrate  towards  the  electrode  with  the  opposite  charge. 
Its  velocity  however,  is  slowed  down  by  the  drag  action  of  the 
surrounding  double  layer,  which  is  pulled  by  the  field  in  the 
opposite  direction. 


The  mobility  of  the  particle  in  an  unbounded  fluid  is  derived  by 
equating  the  electric  force  with  the  frictional  resistance  and  the 
relaxation  force.  Consequently,  the  mobility  u  is  given  bv: 

u  =  Ee2/4h1)  [l-rf(kr)] 

.  f  U 

where : 

u  -  mobility 
E  -  field  strength 
£  -  dielectric  constant  of  fluid 
J?  -  Zeta  potential  of  particles 
-  viscosity  of  fluid 
k  -  1/double  layer  thickness 
r  -  radius  of  particle 

f(kr) .which  varies  between  0  and  1  for  small  and  large  values  of  kr 
respectively,  represents  the  relaxation  phenomenon,  for  k<<<r,  the 
above  expression  reduces  to  the  Smoluchowski  equation,  namely: 

u  =  (E  x  £  x  1  )/ 

u  varies  between  0  and  20xl04  cm/sec/lv/cm. 

The  electorphoretic  yield  is  given  by: 

f 

y  =  a  |uECSdt 
J 

where  a  -  yield  constant 

C  -  cone,  of  particles. 

S  -  area  of  electrode 
t  -  time 

In  the  above  equations  it  is  assumed  that  the  contribution  of  the 
hydrodynamic  ”elocity  to  deposition  is  negligible  its  role  being 
confined  to  maintenance  of  the  suspension. 

As  for  the  phenomena  that  take  place  at  the  electrode,  once  the 
particles  reach  it,  three  possibilities  are  considered. 

One  hypothesis  is  that  the  particles  undergo  an  electrode  reaction 
which  neutralises  them;  this,  however,  does  not  account  for  the 
fact  that  deposits  were  obtained  on  a  polymeric  diaphraqm  located 
between  two  electrodes  f4). 


The  second  hypothesis  states  that  the  particles  are  brought  to  the 
electrode  by  a  field  which  exerts  sufficient  force  to  overcome 
their  mutual  repulsion,  thus  allowing  them  to  come  close  enough  for 
the  London-Van  der  Waals  forces  of  attraction  to  predominate  f5J. 
According  to  the  Verwey  &  Overbeek  [5],  the  minimal  field  strength 
necessary  for  this  purpose  (calculated  from  the  energy  of  particles 
interaction)  is: 

E  =  Z?/3zj  r 

A  third  hypothesis  asssumes  that  secondary  processes  taking  place 
at  the  electrode  can  produce  ions  which  coaqulate  the  particles  by 
discharging  them,  or  produce  hydroxides  which  polymerize  and  adsorb 
on  particles,  thus  holding  them  together  (6). 

In  the  previous  discussion  the  assumption  was  that  the  particle  is 
suspended  in  an  unbounded  fluid.  When  electrophoresis  in  porous 
structures  is  considered,  the  presence  of  -rigid  boundaries  affects 
both  the  electric  field  and  the  velocity  of  the  particles. 

According  to  J.  Anderson( 7 ,  8,  9)  the  effect  of  the  pore  wall  on 

particle  velocity  is  threefold: 

-  The  applied  electric  field  exerts  a  force  on  the  d.l.  at  the 

pore  wall  if  the  latter  is  charged;  elect! oosmotic  flow  or  the 
fluid  is  thus  produced,  which  either  augments  or  opposes  the 
electrophoretic  velocity  of  the  particle,  depending  on  the 
polarity  of  ^  vs ^  . 

-  The  pore  wall  distorts  the  electric  flux  around  the  particle, 
thereby  intensifying  the  local  electric  field  so  that  particle 
velocity  is  augmented. 

-  The  pore  wall  creates  additional  viscous  stresses  in  the  fluid 
which  slow  down  the  particle. 

A  more  detailed  analysis  of  electrophoretic  penetration  is  given  in 
another  report  on  this  study  (10). 
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Parameters  vMch  affect  electrophoretic  deposition  and  penetration 
are : 

-  Zeta  potential  of  the  particles 

-  Particle  size 

-  Particle  concentration 

-  Dielectric  constant  and  viscosity  of  solvent 

-  Electric  field  strength. 

-  Time 

In  this  work  the  feasibility  of  electrophoretic  deposition  of 
several  kinds  of  ceramic  particles  was  first  evaluated,  after  which 
the  effect  of  several  deposition  parameters  on  the  amount  of 
colloidal  Si02  induced  into  the  pores  of  the  substrate  was 
studied. 


Experimental  Procedure 


The  deposition  and  impregnation  experiments  were  carried  out  on 
porous  graphite  substrates. 

The  graphite  was  UCAR  Grade  45,  'with  481  volume  porosity  and  an 
average  pore  size  of  60u  .  Specimens  of  20  x20  x  7mm  were  mostly 
used.  The  ceramic  powders  deposited  were: 

-  colloidal  Si02  -  Pyrogenic  Aerosil,  submicron  with  a  specific 

-  area  of  450  m2/g 

-  fused  Si02  -  l-40u 

-  SiC  -  <40ii 

-  SiN  -  0.2-1. Op 

The  ceramic  particles  were  suspended  m  water,  isopropanol,  ethanol 
and  pentanol. 


The  particle  concentrations  were  3,  13  and  30q/l. 
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Deposition  was  cairied  out  at  constant  voltage  with  the  current 
decreasing  with  time.  Electric  field  strength  was  varied  between  5 
and  150  v/cm  and  deposition  time  between  5  and  120min. 

Prior  to  deposition,  specimens  were  rinsed  with  acetone  for  5  min 
in  an  ultrasonic  bath.  After  deposition  they  were  dried  in  air  for 
24  hrs  before  examination.  For  quantitative  determination  of  the 
impregnated  material  (colloidal  Si02  only)  specimens  were  weighed 
before  and  after  deposition  to  an  accuracy  of  O.lmg.  In  this  case 
the  external  deposit  was  removed  physically  before  weighing.  As  a 
cross  check  the  amount  of  Si02  induced  into  the  porous  structure 
was  determined  by  weighing  the  residue  after  burn-off  of  the 
graphite  at  950°C  for  20hrs. 

The  specimens  were  examined  by  optical  and  electron  microscopy  both 
on  the  surface  and  in  crossection.  The  crossections  were  obtained 
by  fracturing  the  specimen  after  removal  of  the  external  deposit. 


RESULTS 

A  macroscopic  view  of  deposits  of  Sic,  SiN  and  fused  Si02  is  seen 
in  Fiq.  2.  The  deposits  are  in  a  ’’green”  state ( unsintered ) 
and  since  no  special  precautions  had  been  taken,  the  deposits 
cracked  during  dryinq.  The  colloidal  Si02  formed  a  transparent 
deposit  and  is  therefore  presented  at  larger  magnification  in  the 
SEM  (fig  3  ).  which  shows  a  fluffy  feather-like  substance. 

The  presence  of  colloidal  Si02  in  the  pores  of  the  graphite  was 

demonstrated  by  examining  the  crossection  of  a  specimen  by  SEM 
(fig. 4)  and  mapping  the  Si  (fig  5)  on  the  same  area.  Clusters  of 
colloidal  Si02  are  seen  inside  the  crossection,  the  distribution 
of  elemental  Si  coinciding  with  the  clusters. 

Fig  6  shows  the  skeleton  of  Si02  which  remained  after  burn-off  of 

the  graphite  at  900°C  for  2  hrs.  Presence  of  the  Si02  in  the 

crossection  of  the  skeleton  is  seen  m  fiq.  7. 


V 
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Quantitative  data  of  the  penetrated  Si02  as  function  of  solvent 
properties  and  deposition  parameters,  are  given  m  Figs.  8-11.  The 
amount  of  Si02  is  expressed  as  weight  percent  of  the  specimen  . 

The  amount  of  penetrated  Si02  as  function  of  the  ratio  between 
dielectric  constant  and  viscosity  of  the  solvent  is  shown  m  Fig.  8 
for  constant  field,  Si02  concentration  m  solvent  and  deposition 
time.  A  strong  influence  of  the  solvent  properties  is  revealed. 

The  dependence  on  the  field  intensity  for  three  different  Si02 
concentrations  in  water  and  for  two  concentrations  in  propanol  is 
given  in  Fig.  9.  To  each  concentration,  there  corresponds  an 
optimal  field.  The  Si02  concentration  m  the  solvent  affects  the 
amount  of  penetrated  Si02  for  otherwise  identical  conditions; 
this  is  seen  in  Fig.  10  for  deposition  from  water  and  from 
propanol . 

The  effect  of  deposition  time  on  the  amount  of  penetrated  Si02  is 
shown  in  Fig.  11  again  for  deposition  from  water  and  from  propanol. 

Discussion 

The  feasibility  of  depositing  oxide  and  non-oxide  ceramics  on  an 
electrically  conductive  substrate  by  electrophoresis  is 
demonstrated  m  Fig.l. 

Charging  of  the  Si02  particles  is  probably  the  consequence  of 
dissociation  of  silanol  groups  in  the  silica: 

-  SiOH  +  OH  - >=  SiO-  +  H20 

The  charge,  thus  obtained,  is  negative  and  in  fact  the  particles 
migrate  towards  tne  anode  in  the  electric  circuit.  Charging  of 
non-oxide  ceramics  takes  place  probably  via  an  absorption 
mechanism,  with  SiC  acquiring  a  negative  charge  and  depositing  on 
the  anode  and  SiN  acquiring  a  positive  charge  and  depositing  on  the 
cathode . 
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The  cracks  in  the  coating  formed  during  drying  as  a  result  of 
shrinkage.  A  study  is  currently  under  way  to  modify  the  drying 
procedure,  so  that  contraction  stresses  will  be  minimized  and 
cracking  prevented. 

In  order  to  verify  the  possibility  of  induction  of  ceramic 
particles  into  the  pores,  crossections  of  coated  graphite  were 
prepared  by  fracturing  the  specimens.  The  presence  of  Si02 
clusters  m  the  crossection  is  evidenced  both  morphologically  in 
fig  3.  and  by  microanalysis  in  fig  4.  The  clusters  vary  in  size 
possibly  reflecting  the  variation  of  pore  size.  The  distribution  of 
elemental  Si  coincides  as  expected,  with  cluster  locations. 

A  most  interesting  result  was  obtained  during  quantitative 
determination  of  the  impregnated  Si02.  it  was  found  that  in 
specimens  in  which  the  induced  Si02  was  at  the  higher  levels 
(such  as  0.6-2.51)  a  Si02  skeleton  of  the  same  shape  and 
dimensions  as  the  original  specimen  remained  after  removal  of  the 
graphite . 

Moreover,  the  presence  of  the  Si02  throughout  the  crossection  of 
the  skeleton  indicates  impregnation  of  the  whole  crossection  of  the 
specimen . 

The  effect  of  solvent  properties  on  penetration  can  be  predicted 
from  the  theoretical  analysis  in{10) .  Penetration  is  enhanced  at 
large  Peclet  numbers: 

Pe  =  Ub ,  where  U  is  the  particle  velocity,  b  the  mean  pore 
D 

radius  and  D  the  diffusion  coefficient  of  the  particle.  The 
electrophoretic  velocity  is  proportional  to  e/£  of  the  solvent. 

Four  solvents  were  tested  with  e/ij^  ratios  ranqing  from  4.2  for 
pentanol  to  81  for  water.  As  predicted,  the  amount  of  penetrated 
Si02  increased  with  this  ratio  (fig  /). 
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The  electric  field  has  a  dual  effect  on  the  extent  of  particle 
penetration.  Due  to  increase  of  the  particle  velocity  and  Peclet 
number  with  field  strength,  penetration  is  enhanced  by  it.  However, 
the  field  enhances  also  the  build-up  of  an  external  deposit  on  the 
substrate  which  blocks  penetration.  Therefore,  there  exists  an 
optimal  field  strength  (fig.  8),  its  value  decreasing  with 
increased  particle  concentration,  because  of  faster  build-up  of  the 
external  coating  at  higher  concentrati ons . 

Penetration  increases  with  particle  concentration  both  in  water  and 
propanol.  At  low  field  strengths  (5V/cm)  the  relationship  is  linear 
for  the  range  of  concentrations  studied  (fig. 9). 

Penetration  increases  with  deposition  time  but  reaches  a  plateau 
after  a  period  which  depends  on  solvent  concentration  and  field 
strength.  The  fact  that  the  plateau  is  reached  before  full 
impregnation  of  the  pores,  is  probably  due  to  the  blockage  by  the 
external  coating.  Thus  the  extent  of  penetration  could  be  increased 
if  a  way  is  found  to  prevent  surface  deposition. 

Conclusions 

A  variety  of  ceramic  materials,  both  oxides  and  nonoxides,  were 
deposited  electrophoretically  on  a  porous  graphite  substrate. 

In  addition  to  formation  of  surface  deposits,  electrophoretic 
induction  of  ceramic  particles  into  the  porous  substrate  was 
demonstrated. 
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Quantitative  studies  of  the  amount  of  penetrated  material  as 
function  of  deposition  parameters,  were  carried  out  for  colloidal 
Si02.  it  was  shown  that  penetration  is  enhanced  by  a  high 
dielectric  constant  and  low  viscosity  of  the  solvent,  by  particle 
concentration  and  by  field  strength.  Due  to  simultaneous  build-up 
of  an  external  deposit,  a  plateau  is  reached  in  the  extent  of 
penetration  as  function  of  time,  prior  to  filling  of  the  whole  pore 
volume.  Prevention  or  inhibition  of  surface  deposition  would  favor 
an  increased  extent  of  penetration. 

However,  at  the  optimal  conditions  reached  so  far,  penetration  of 
the  whole  crossection  of  the  specimens  was  demonstrated. 
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Fig.l  Egg-shaped  ionic  cloud  and  forces  acting  on  a  moving 
particle 


Fiq.2  "Green-'  deposits  of  SiC,  SiN  and  fused  Si02  (from  left  to 
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Fig. 3  Colloidal  SiOa  deposit  on  graphite 
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Fig. 5  X  -  ray  mapping  for  Si  on  area  seen  in  Fig.  4 
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Fig. 8  Amount  of  Si02  that  penetrated  the  porous  graphite  irom 
various  solvents. 
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Fig. 11  Effect  of  deposition  time  on  amount  of  penetrated  Si02 
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The  electrophoretic  penetration  of  colloidal  silica  particles  into  a  porous  graphite 
substrate  is  investigated.  The  graphite  substrate  is  immersed  in  a  suspension  containing  the 
particles  and  positioned  between  two  electrodes.  An  electric  potential  gradient  between  the 
graphite  and  the  electrodes  is  used  to  drive  the  silica  particles  into  the  pores.  Three  driving 
mechanisms  are  identified:  the  hydrodynamic  drag  force  exerted  on  the  particles  due  to  the 
electroosmotic  flow  of  the  solvent  inside  the  pores,  the  electrophoretic  force  exerted  on  the 
particles  and  the  stochastic  Brownian  force  due  to  thermal  fluctuations  of  the  solvent 
molecules.  While  subjected  to  these  forces,  the  particles  may  reach  the  walls  of  the  pore 
and  the  short  range  van  der  Waals  forces  may  cause  their  capturing  and  deposition  onto  the 
walls. 


The  objectives  of  this  paper  are:  a)  to  predict  the  penetration  depth  of  a  single 
ceramic  particle  moving  inside  a  porous  substrate  under  the  effect  of  an  electric  potential 
gradient,  b)  to  derive  the  nondimensional  parameters  characterizing  the  motion  of  the 
ceramic  particles, and  c)  to  gain  a  physical  insight  on  the  various  mechanisms  governing 
penetration.  Qualitatively,  the  results  are  that  penetration  depths  are  governed  by  a 
favorable  (if  large)  Peclet  number  and  unfavorable  (if  large)  local  deposition  rate  number. 
Quantitative  results  are  also  provided. 
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Deep  penetration  and  deposition  of  ceramic  materials  on  oxidizable  porous  substrates 
has  been  a  major  concern  of  the  aero-space  industry.  The  extremely  strong,  light  and 
sometimes  porous  composite  of  carbon-carbon  may  undergo  undesired  processes  if 
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exposed  to  a  high  temperature  oxidizing  environment.  Chemical  Vapor  Deposition  (CVD) 
of  itett  silica,  silicon-carbide  etc.  onto  the  walls  of  the  pores,  has  been  suggested  in  the 
past  (1-2)  to  protect  the  composite  from  deterioration.  Alas,  CVD  has  proven  to  be 
extremely  slow  and  only  shallow  penetration  has  been  detected.  A  rather  new  method  to 
protect  porous  substrates  by  electrophoresis  is  studied  and  described  in  detail  in  part  H  In 
essence,  an  electric  potential  gradient  is  used  to  drive  colloidal  particles  deep  into  the  voids 
of  a  porous  model  material  for  carbon-carbon  -  porous  graphite.  The  purpose  of  this  paper 
(Part  I)  is:  a)  to  provide  an  overview  of  the  existing  theory  and  mathematical  models 
closely  related  to  Deep  Electrophoretic  Deposition  (DED),  b)  to  suggest  a  theory  and  a 
mathematical  model  by  which  DED  can  be  analyzed  and  the  non-dimensional  parameters 
governing  the  process  be  obtained,  and  c)  to  predict  penetration  depths  for  a  given  set  of 
parameters. 

The  motion  of  a  single  particle  under  the  effect  of  electrophoretic  forces  has  first 
been  addressed  by  Smoluchowski  (1918)  and  subsequently  by  many  others  (3-15). 
Smoluchowski  predicted  that  a  rigid  spherical  particle  possessing  an  electric  double  layer 
and  embedded  in  an  unbounded  flow  field  would  be  forced  to  move  if  subjected  to  an 
electrical  potential  gradient  The  mobility  of  the  particle  depends  linearly  upon  the  dielectric 
constant  of  the  fluid,  the  potential  gradient,  the  zeta  potential  of  the  particle  and  is  inversely 
proportional  to  the  fluid  viscosity.  Wall  effects  were  accounted  for  by  Morrison  &  Stukel 
(1970)  and  Keh  &  Anderson  (1985)  for  the  case  of  a  particle  travelling  in  close  proximity 
to  the  containing  boundaries  of  the  flow  field.  Keh  &  Anderson  (1985)  have  shown  that  (if 
electroosmosis  is  disregarded)  the  particle  would  experience  an  increasing  drag  force  and 
reduced  mobility.  It  was  also  shown,  that  due  to  the  electric  field,  the  insulating  rigid  walls 
would  induce  an  electroosmotic  flow,  its  direction  depending  on  the  sign  of  the  wall  zeta 
potential.  Thus,  a  neutral  particle  immersed  in  the  fluid  would  be  dragged  by  the 
electroosmotic  flow  and  move  with  almost  identical  velocity.  In  case  the  particle  is  also 
subjected  to  an  electrophoretic  force,  the  particle  would  move  with  the  combined 
electrophoretic  and  electroosmotic  velocities.  Superimposition  is  allowed  due  to  the 
linearity  of  the  governing  low  Reynolds  number  flow  field  equations.  In  addition,  a  sub- 
micron  particle  would  experience  an  erratic  Brownian  motion  due  to  the  thermal 
fluctuations  of  the  fluid  molecules.  Consequently,  deterministic  evaluation  of  the  particle 
path  under  the  combined  effects  of  electroosmosis,  electrophoresis  and  Brownian  motion  is 
invalid  and  a  probabalistic  approach  must  be  adopted  (van-Kampen  (1983)). 
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The  fluid  motion  inside  a  porous  substrate  under  the  effect  of  a  pressure  gradient 
(known  as  Darcy's  law)  has  been  addressed  by  many  investigators.  To  devise  the  micro- 
pattern  of  the  flow  field,  unit  cell  models  (19-21)  as  well  as  more  advanced  finite  elements 
approaches  for  spatially  periodic  structures  (22-27)  were  applied.  However,  no  equivalent 
treatment  exists  for  the  case  of  electroosmodc  flow  through  a  porous  structure,  namely,  the 
flow  generated  under  the  effect  of  an  electric  potential  gradient  Since,  a  particle  would  be 
dragged  by  this  interstitial  electroosmotic  flow,  its  evaluation  is  of  foremost  significance  if 
one  desires  to  calculate  mean  penetration  depths  of  the  ceramic  particles.  A  general 
approach  in  which  we  modify  Darcy's  law  so  that  it  would  incorporate  electroosmotic 
effects  is  addressed  in  a  separate  paper  (Haber  (1990)). 

The  second  effect  to  be  accounted  for  is  the  direct  electrophoretic  force  exerted  on  a 
particle  travelling  inside  the  small  pores.  Numerous  treatments  exist  for  the  motion  of  a 
small  rigid  particle  at  low  Reynolds  numbers,  wall  effects  excluded  (29)  and  included  (29- 
32).  The  fine  structure  of  a  porous  substrate  enveloping  a  particle  was  treated  in  various 
ways,  for  instance,  a  spherical  particle  travelling  along  an  infinitely  long  cylinder  of  circular 
ctoss-section,  a  spherical  particle  embedded  in  a  flow  field  which  is  bounded  by  a 
spherical  envelope  with  pre-assumed  free  surface  (or  free  vorticity)  boundary  conditions, 
and  a  small  spherical  particle  moving  inside  a  spatially-periodic  lattice  of  pre-arranged  large 
spherical  particles.  Notwithstanding,  all  of  the  previous  models  assumed  no-slip  boundary 
conditions  over  the  particle  surface,  whereas  the  electrophoretic  force  is  applied  via  a 
unique  slip  boundary  condition.  Thus,  the  known  approaches  must  be  rectified 
accordingly. 

* 

The  third  effect,  the  random  Brownian  force  due  to  the  thermal  fluctuations  of  the 
fluid,  was  treated  by  two  different  methods  in  the  past.  The  first  method  applies  a 
Lagrangian  viewpoint  of  the  problem  utilizing  Langevin's  equation  to  evaluate  mean 
particle  velocity  and  dispersion  (33,34).  The  second  method  applies  an  Eulerian  approach 
where  a  Fokker-Planck  equation  is  formulated  and  a  moment  method  developed  by  Taylor- 
Aris  is  utilized  to  derive  the  pertinent  means  (35-37).  Generally,  it  is  accepted  (and  for 
linear  cases  can  be  proven  (18))  that  these  methods  would  yield  identical  results  and  the 
choice  between  the  methods  is  a  matter  of  convenience. 

In  the  next  chapter  we  shall  address  the  three  foregoing  mechanisms  where  we 
provide:  a)  a  short  recapitulation  of  the  basic  equations  governing  electroosmotic  flow  in 
porous  structures  derived  elsewhere  (Haber  1990),  b)  a  model  for  the  motion  of  a  particle 
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in  a  long  cylindrical  tube  under  the  effect  of  an  electric  field,  (a  model  simulating  the  motion 
of  a  particle  in  a  porous  structure),  and  c)  a  model  for  the  stochastic  behavior  of  the 
Brownian  particle  under  the  combined  effect  of  electroosmosis  and  electrophoresis  utilizing 
Fokker-Planck's  equation. 
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,1a. Electroosmosis  in  Porous  Media 

The  flow  through  porous  media  generated  by  an  electric  potential  gradient  was 
analytically  investigated  by  Kozak  &Davis  (1986,1989)  using  a  unit  cell  model  and  by 
Haber  (1990)  using  a  simple  approach  based  on  the  known  solution  of  the  flow  field 
through  a  long  circular  tube  (in  close  similarity  to  one  of  the  approaches  used  to  analytically 
prove  Darcy's  law).  Similar  approaches  for  pressure  driven  flows  were  used  to  determine 
filtration  efficiencies  in  porous  substrates  (39).  We  recapitulate  here  the  basic  assumptions 
and  results  presented  by  Haber  (1990). 

The  governing  field  equations  for  the  flow  through  a  long  circular  tube  induced  by 
an  electric  potential  gradient  are  based  on  Stokes'  equations  for  low  Reynolds  number 
flows,  i.e. 

7]V2v  =  Vp  V-v  =  0  [1] 

where  v  and  p  are  the  velocity  and  pressure  fields  and  V  stands  for  tire  viscosity  of  the 
fluid.  The  boundary  conditions  which  the  velocity  field  satisfies  over  the  cylinder  walls  are 
(17): 


‘4  Trqdz 


v  =0 


@  =  b 


m 


where  e  is  the  dielectric  constant  of  the  fluid,  is  the  zeta  potential  at  the  cylinder  wall,  <p 
stands  for  the  electric  potential,  b  is  the  radius  of  the  tube  and  (r,z)  are  polar  coordinates, 
the  z  coordinate  coinciding  witli  the  cylinder  axis.  It  should  be  noted  that  die  no-slip 
condition  can  no  longer  be  applied  and  in  its  place  we  use  condition  [2] . 
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The  solution  for  v  is  easily  derived  for  a  fully  developed  flow,  namely,  dvjdz=0. 
It  is  a  simple  superposition  of  the  effects  of  the  pressure  gradient  and  the  electric  potential 
gradient.  Integrating  over  the  tube  cross-section,  the  mean  interstitial  velocity  is  obtained: 

v  .  p] 

‘  St?  dz  47tr?  dz 

Integrating  over  the  entire  cross-section  of  the  porous  medium  the  superficial 
velocity  is, 

v  _  ep<b</b2>dp  e,eCwdtj>  [4] 

'  87?  dz  4777?  dz  '  J 


where  ep  stands  for  the  medium  porosity. 

Since  a  homogeneous  porous  medium  contains  pores  in  all  possible  directions, a 
similar  procedure  can  be  applied  for  pores  oriented  along  the  x  and  y  directions. 
Expressions  similar  to  [4]  can  be  obtained  for  v*  and  vy  and  further  generalized  to  derive 
the  following  expression  for  v  in  vectorial  notation, 

v  =  -rV/7-y70  [5] 

where 


[6] 


Here  F  and  yean  be  viewed  as  the  hydrodynamic  and  electric  permeabilities, 
respectively.  Equation  [5]  can  be  viewed  as  a  generalization  of  Darcy's  law  for  the  case 
where  an  electrical  potential  gradient  as  well  as  a  pressure  gradient  is  imposed  over  a 
porous  medium.  It  is  obvious  that  these  two  forces  are  decoupted  and  a  potential  gradient 
alone  can  induce  electroosmotic  flow  inside  a  porous  structure.  Increasing  the  electrical 
permeability,  y,  would  cause  an  increase  in  the  electroosmotic  flow.  Thus,  an  increase  in 
medium  porosity,  fluid  dielectric  constant  and  zeta-  potential  of  the  pore  walls  and  decrease 
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in  fluid  viscosity  would  cause  enhanced  flow  through  a  porous  structure  for  a  given 
potential  gradient. 

In  addition,  one  must  satisfy  mass  conservation  which  results  in  the  following 
continuity  equation  for  incompressible  fluid,- 

V- v  =  0.  [7] 

Equations  [5]  and  [7]  and  the  appropriate  boundary  conditions  on  <p  and  p  do  not 
constitute  a  complete  set  by  which  the  velocity,  pressure  and  potential  fields  can  be 
evaluated.  An  additional  relation  between  the  flux  of  the  electric  charge  j  and  the  electric 
potential  and  pressure  gradient  can  be  established  in  a  manner  similar  to  the  derivation  of 
equation  [5] : 


j  =  -yVp-KV<t>  [8] 

where  K  is  the  electric  conductivity.  It  is  interesting  to  note  that  y  appears  both  in  [5]  and 
[8]  in  full  agreement  with  Onsager's  reciprocal  rule.  The  exact  evaluation  of  K  is  of  no 
interest  at  the  moment  and  will  not  be  addressed.  Similat  to  [7]  we  obtain, 

V-j  =  0  [9] 

which  is  simply  Kirchoffs  law  applied  to  a  porous  medium. 

Utilizing  [5],  [7],  [8]  and  [9]  it  is  easy  to  show  that  <j>  and  p  can  be  decoupled  and 
satisfy  the  following  equations  for  a  homogeneous  medium: 

V20  =  V2p  =  0  ■  [10] 

i.e.  <p  and  p  are  harmonic  functions  which  can  easily  be  derived  for  a  given  set  of 
boundary  conditions. 

For  the  very  simple  case  of  an  infinite  slab  of  thickness  c,  which  is  subjected  to  a 
potential  drop  #  and  a  zero  pressure  gradient  (along,  say,  the  z  direction),  the  steady- 
state  solutions  for  <p  and  v  are  simply 
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<P  =  <Pi+(<p2-<Pi)z/c  ,  v=u*=- r(4>i-Mc  ■  [ii] 

Thus,  the  electroosmotic  velocity  through  the  slab  is  uniform  and  constant.  We 
adopt  the  symbol  U*  for  the  velocity  solution  of  this  simple  case  and  shall  apply  it  later  on 
in  chapter  3c. 

lii.  Electrophoretic.  Motion  of  a  Particle 


The  electrophoretic  velocity  U"  of  a  small  rigid  spherical  particle  suspended  in  an 
unbounded  flow  field  of  viscosity  7?  and  subjected  to  an  external  electric  potential  gradient 
E  =  -70  was  first  obtained  by  Smoluchowski, 

'121 

where  (p  is  the  zeta  potential  over  the  particle  interface.  The  basic  assumptions  made  to 
derive  [121  are  that  the  Debye  screening  length  is  much  smaller  than  the  particle  dimension 
and  that  the  particle  translates  in  an  unbounded  flow  field. 

If  the  particle  translates  at  close  proximity  to  a  rigid  surface,  wall  effects  must  be 
accounted  for  (e.g.  Keh  and  Anderson(1985)).  Notwithstanding,  for  particle  dimensions 
which  are  more  than  ten  times  smaller  than  its  distance  to  the  wall,  Smoluchowski’s 
equation  is  an  excellent  approximation  (given  that  no  electroosmotic  flow  exists). 

It  is  a  well  known  result  from  the  theory  of  low  Reynolds  number  flows  (40,41) 
that  a  particle  approaching  a  rigid  wall  (subjected,  for  instance,  to  gravity  forces)  will 
experience  a  growing  hydrodynamic  resistance  which  at  the  limit  of  zero  gap  increases  to 
infinity.  However,  the  condition  of  slip  velocity  over  the  particle  boundaries  will  cause  that 
resistance  to  grow  in  a  slower  manner  up  till  a  gap  of  the  order  of  the  Debye  screening 
length.  An  exact  treatment  of  this  hydroaynamically  singular  behavior  awaits  exploration. 
However,  over  these  very  small  distances  one  has  to  consider  other  surface  (e.g.  Ixmdon- 
van  der  Waals)  forces  which  eventually  lead  to  the  attraction  of  the  particle  to  the  walls. 
Thus,  it  seems  reasonable  to  assume  that  Smoluchowski’s  equation  can  be  used  as  a  first 
order  approximation  for  the  mobility  of  a  particle  over  the  entire  range  of  particle-wall 
gaps. 
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Henceforth,  we  shall  assume  that  equation  [12]  is  applicable  for  the  case  of  a  single 
rigid  particle  immersed  in  a  fluid  of  viscosity  Tj  and  bounded  by  the  wails  of  a  long  tube. 


3c.  The  Probability  Distribution  for  Particle  Penetration  Depths 


3c.  1  Statement  of  The  Problem 

Small  sub-micron  particles  travelling  inside  the  pores  of  a  porous  substrate  (see 
figure  1)  are  strongly  affected  by  the  thermal  fluctuations  of  the  fluid  molecules  and 
experience  Brownian  motion.  This  stochastic  motion  must  be  superimposed  to  the 
deterministic  electrophoretic  motion  and  electroosmotic  induced  velocity  which  are 
primarily  parallel  to  the  direction  of  the  pore  axis.  The  Brownian  motion  causes  the  particle 
to  cross  streamlines  and  sample  all  possible  transverse  positions  which  leads  to  the 
phenomenon  known  as  Taylor  dispersion. 
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Tracing  down  exactly  the  trajectory  of  a  particle  is  therefore  of  no  consequence.  What  we 
seek  is  the  probability  distribution  that  a  particle  entering  the  porous  substrate  would  reach 
a  given  depth  and  not  be  deposited  on  the  walls  during  its  erratic  motion  inside  the  pores. 

In  addition,  the  very  complex  and  random  structure  of  the  porous  substrate  makes  it 
impossible  to  obtain  more  than  a  formal  representation  of  the  equations  and  boundary 
conditions.  The  mathematical  formulation  must  hinge  on  a  simplified  geometrical  model  for 
the  porous  structure  which  would  be  amenable  to  mathematical  analysis.  One  of  the  very 
common  approaches  is  to  treat  a  single  pore  as  a  long  circular  tube  of  mean  pore  diameter 
its  axis  coplanar  with  the  local  superficial  velocity  direction.  Such  a  model  is  handicapped 
by  the  fact  that  no  pronounced  lateral  dispersion  of  the  particles  across  the  porous  substrate 
is  allowed.  However,  since  only  the  mean  longitudinal  particle  motion  is  sought  such  a 
model  is  expected  to  provide  valid  results  on  the  main  parameters  affecting  penetration 
depths.  Figure  2  describes  the  basic  geometrical  and  kinematical  parameters  applicable  to 
the  problem. 


The  general  differential  equations  governing  tire  problem  is  Fokker-Planck’s 
equation  for  the  probability  distribution.  In  our  case  it  assumes  the  form, 


dp 

dt 


+  V,- J,  +  Vr- Jr  +  Ve- J6 


-8(z)5(r-r0)8(e-80)8(t) 


[13] 
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where 


J,  =  U p-D,Vsp,  J  r  =  ~DrVrp,  Je  =  -DeVep 

D=kTM,  D=kTMr  De  =  kTMe  [14] 

U=U‘(f)  +  U"(f) 

Here  p(r,z,6,t  /  r0,0,60)  stands  for  the  conditional  probability  that  a  panicle  initially 
located  at  (r0,0,60)  reaches  the  infinitesimal  volume  around  ( r,z,G)  after  time  t.  The  polar 
coordinate  system  ( r,z,6)  coincides  with  the  pore  axis  where  r  stands  for  the  radial 
distance  from  the  pore  axis,  z  is  the  distance  measured  along  the  pore  axis  from  the  pore 
entrance  and  6  is  the  azimuthal  angle.  The  symbol  1)  stands  for  the  determinisdc  pan  of 
the  velocity  of  the  particle,  composed  of  its  electrophoretic  velocity  and  the  clectroosmoiic 
velocity  of  the  suspending  field.  It  must  be  stressed  that  the  interstitial  flow  due 
electroosmosis  U"  does  not  possess  the  Commonly  agreed  upon  parabolic  (Poiseuille) 
form.  It  is  very  well  approximated  by  plug  flow  (as  shown  in  figure  2)  for  pore  diameters 
larger  than  Debye's  screening  length  (a  condition  normally  met).  Thus,  U  is  a  uniform 
function  with  respect  to  the  spatial  coordinates  as  long  as  the  particle  diameter  is  at  least  ten 
times  smaller  than  the  pore  diameter.  The  flux  Jz  parallel  to  the  axis  of  the  pore  depend  on 
particle  convection  and  diffusion  along  the  pore  axis,  whereas  tire  fluxes  Jt  and  Je  arise 
from  diffusion  only.  The  diffusion  coefficients  Dv  DT  and  D0,  the  first  one  parallel  and  the 

latter  perpendicular  to  the  pore  axis,  depend  on  the  absolute  temperature  T  of  the  liquid  and 
the  respective  particle  mobilities  Mz  Mx  and  Mq.  This  generalized  Stokes-Einstein 

relationship  provided  in  [14]  accounts  for  wall  effects- through  the  anisotropic  form  of  the 
mobility  tensor.  From  the  physical  point  of  view,  it  is  obvious  that  the  particle  would 
experience  unequal  hydrodynamic  resistance  to  its  motion  perpendicular  and  parallel  to  the 
pore  axis.  The  mobility  tensor  components  were  approximately  evaluated  in  Happel  and 
Brenner  (1973).  It  is  shown  that,  for  panicles  of  radius  a  small  compared  to  their  distance 
to  the  wall,  the  mobilities  Mz,  Mr  and  Mq  are  approximately  equal  and  possess  the 
following  form,  ' 


Mt~  Mr~  Me  ~  \/(67trja)  . 


[15) 
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The  boundary  condition  which  must  be  satisfied  by  the  probability  distribution  p 
is, 

-D^-  =  Kp  @  r  =  b  [16] 

dr 

where  k  is  the  local  deposition  rate  and  b  stands  for  the  pore  mean  radius. 

The  left  hand  side  of  [16]  represents  the  diffusional  flux  of  probability  towards 
the  pore  walls  whereas  the  right  hand  side  represents  the  probability  that  a  particle  located  at 
close  proximity  to  the  walls  would  indeed  be  deposited.  The  numerical  value  of  K  depends 
on  particle  shape  and  dimension,  electrostatic  and  hydrodynamic  forces  and  the  physical 
properties  of  the  materials  comprising  the  particles,  the  porous  substrate  and  the 
suspending  liquid.  Indeed,  it  is  extremely  difficult  to  obtain  k  from  first  principles.  A 
feasible  method  to  attain  this  end  would  be  to  devise  an  experiment  in  which  one  would 
utilize  the  given  materials  of  the  ceramic  particles,  porous  substrate  and  liquid. 


It  is  very  useful  to  state  the  non-dimensional  parameters  and  time  scales  of  the 
problem  before  we  attempt  to  obtain  an  analytical  or  numerical  solution  of  [13].  It  is  easy  to 
show  from  [6]  that  the  pertinent  time  scale  within  which  steady  electroosmotic  flow  can  be 
established  is, 


[17] 


where  E  is  the  externally  applied  electrical  field. 

Similarly,  it  can  be  shown  from  [13],  [14]  and  [16]  that  three  time  scales  exist  for 
the  penetration  process  of  Brownian  ceramic  particles  into  the  porous  substrate.  The  first 
time  scale  is  linked  with  the  convective  term  of  the  Fokkcr-Planck  equation,  namely, 


c 


[18] 
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which  determines  the  time  it  would  take  the  particle  to  cross  the  porous  slab. 

The  second  time  scale  is  related  to  the  diffusive  process,  namely,  the  time  it  lakes  a 
Brownian  particle  to  cross  streamlines  and  sample  all  radial  positions  inside  the  tube, 


The  third  time  scale  is  obtained  from  boundary  condition  [16]  and  determines  the 
deposition  rate,  namely, 


If  we  assume  that  the  Brownian  particles  were  introduced  into  the  solution  after 
electnoosmotic  steady  state  flow  has  been  established,  (U  is  time  independent)  non- 
dimensionalization  of  the  Fokker-Planck  equation  and  the  associated  boundary  conditions 
results  in  two  non-dimensional  parameters  which  determine  the  solution  completely.  The 
first  parameter  can  be  identified  as  the  Peclet  number  of  the  problem, 

[21] 

D  D 

which  determines  the  relative  significance  of  the  convective  process  vis-a-vis  the  diffusive 
process. 

The  second  nondimensional  parameter,  henceforth  called  as  the  deposition  rate 
parameter,  determines  the  ratio  between  the  diffusion  and  the  deposition  time  scales, 
namely, 


[22] 
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Thus,  large  values  of  X  mean  that  the  deposition  of  particles  on  the  cylinder  walls 
is  governed  mainly  by  diffusion,  whereas  small  values  of  X  mean  that  deposition  is 
controlled  by  the  local  deposition  rate. 


An  analytical  solution  of  [13]  (in  its  non-dimensional  form)  can  readily  be  obtained 
for  the  case  that  U  is  independent  of  the  radial  coordinate  r.  Since  U'  and  U"  are  indeed 
void  of  radial  coordinate  dependence  for  small  particle  radius  to  pore  diameter  ratios,  we 
obtain  after  tedious  manipulations, 


P{riX^)  = 


<r  Z(Jt?'4(<V7V77) 

71=1  Tj’ 


ct-Jckajl)  -alz 

a2+«X(«,.) 


[23] 


[24] 


and/0  is  the  Bessel  function  of  zero  order. 

The  parameters  are  determined  by  the  following  indicial  equation, 


ccJ,(a)  =  Xf0(a)  [25] 

which  has  a  countable  number  of  roots  a  (n=l,2,...). 

n 

The  expression  for  the  probability  distribution  P  in  equation  [23]  is  independent  of 
the  initial  azimuthal  angle  60  and  radial  position  r0 .  This  is  achieved  by  assuming  that  all 
initial  positions  of  the  particle  at  the  pore  entrance  are  equally  probable  and  integrating  the 
conditional  probability  distribution  p  over  all  these  positions.  Thus,  [23]  does  not 
represent  a  solution  for  [13]  it  is  rather  the  analytical  solution  for  the  physically  plausible 
case  based  on  p, 

P(V,C’*)  =  b3jjp(Ti>C,0,  eo)P(Vo<  6o)dA 


[26] 
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2 

where  A0=jtR0  is  the  cross-section  area  of  the  pore  and  P(7)0, 60)  =  1//^  is  the  uniform 
probability  distribution  of  finding  a  particle  at  time  t=0  at  any  arbitraiy  location  at  the  pore 
entrance.  It  is  interesting  to  note  that  P  is  independent  of  6  due  to  the  inherent  axial 
symmetry  of  the  problem. 


Equations  [23]  and  [25]  can  now  be  utilized  to  numerically  compute  more 
physically  meaningful  probabilities.  For  example,  if  we  are  to  inquire  what  is  the 
probability  per  unit  length  that  a  particle  have  deposited  over  the  pore  walls,  a  very  long 
time  after  it  have  been  introduced  into  the  system,  we  simply  calculate, 


dP 

dtj 


dt 


[27] 


where  the  right  hand  side  equality  stems  from  boundary  condition  [16].  It  is  easy  to  show 
that 


P^;F<,A)  =  42£(jr]70(a„r?Vr?') 


a; 


n=l 


(A2 + aX(a.W?+4«1V?+4«?  -  p) 


[28] 


Figure  3  illustrates  the  numerically  computed  values  of  Pd  as  a  function  of  z  for 
Peclet  numbers  1  and  10  and  A  values  0.1, 0.5  and  1.  The  trend  is  very  clear.  High  values 
of  X  cause  the  probability  density  to  decay  very  fast  with  z .  This  means  that  most  of  the 
particles  would  deposit  on  the  walls  very  close  to  the  pore  entrance  and  deep  penetration 
cannot  be  achieved.  Figure  4  illustrates  Pd  as  a  function  of  Peclet  number  Pe  for  X  values 

0.1, 0.5  and  1  and  for  given  pore  depths  (z=0  and  1).  It  manifests  the  very  important  role 
of  particle  convection.  The  probability  density  drops  significantly  with  the  Peclet  number 
Pe  thus  a  more  even  distribution  of  particles  is  expected  along  the  pore  axis  and  deeper 
penetration  is  expected.  In  figure  5  the  probability  density  is  directly  plotted  vs.  the  local 
deposition  parameter  A  for  Peclet  values  0.1, 1.0  and  10  and  two  values  of  z,  showing 
again  the  fast  accumulation  of  particles  near  the  pore  entrance  for  small  values  of  Peclet 
number. 
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Figures  6  to  9  illustrate  the  foregoing  conclusions  vividly.  Based  on  the  probability 
density  Pd  we  can  define  a  new  cummulative  probability 

•It 

P'=  \PM  [29  j 

o 


which  can  physically  be  intetpreted  as  the  long  time  ratio  between  the  number  of  particles 
deposited  over  the  pore  walls  up  till  pore  depth  z/b  and  the  total  number  of  particles 
introduced  at  z=0  and  t=0.  The  figures  elucidate  the  fact  that  deep  penetration  (more  than 
hundred  pore  diameters)  can  be  achieved  only  for  very  low  deposition  parameters  X  ( X  < 
0.1 )  and  for  very  high  Peclet  numbers  (Pe>  100).  Cases  in  which  X  =0  would  obviously 
result  in  the  deepest  penetration.  It  might  be  achieved  by  introducing  a  repelling  force 
between  the  particles  and  the  pore  walls,  stronger  than  van  der  Waals  forces.  The 
asymptotic  value  of  the  cummulative  probability  for  large  values  of  z  is  always  unity,  since 
after  a  long  time  the  probability  that  a  particle  has  been  deposited  anywhere  over  the  pore 
walls  is  unity  (the  particle  cannot  disappear). 

Figure  10  summarizes  the  foregoing  results.  We  define  the  penetration  depth  as  the 
zlb  value  for  which  Pc=0.9 .  In  other  words,  10  percent  of  the  particles  would  penetrate 
deeper  than  the  value  shown  in  the  figure  for  a  given  set  of  parameters,  the  local  deposition 
rate  X  and  the  Peclet  number  Pe.  Obviously,  the  0.9  numerical  value  chosen  for  Pc  is  quite 
arbitrary  and  a  lower  or  a  higher  percentage  could  have  been  chosen.  The  basic 
interpretation  of  the  results,  however,  would  remain  unaltered.  It  is  obvious  from  figure  10 
that  the  penetration  depth  for  a  given  value  of  X  grows  almost  linearly  with  the  Peclet 
number  (which  is  not  at  all  surprising).  The  slope  of  the  curves  depends  strongly  on  X :  the 
smaller  the  X  value  the  deeper  the  particle  penetration.. 


4.  CONCLUSIONS 

a)  Deep  penetration  and  coating  is  feasible  by  electrophoretic  processes. 

b)  Enhanced  penetration  is  obtained  for  large  electrophoretic  Peclet  numbers 
and  small  deposition  rate  constants. 
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c)  Electroosmosis  increases  Peclet  number  and  thus  penetration. 

d)  The  diffusion  coefficient  of  the  colloidal  particles  has  a  dual  effect: 

1)  decreasing  diffusion  coefficient  causes  increasing  Peclet  number  which 
has  a  favorable  effect  on  penetration  depth. 

2)  decreasing  diffusion  coefficient  causes  increasing  deposition  rate  which 
has  an  unfavorable  effect  on  penetration  depth. 

e)  The  process  can  be  controlled  by  a  judicious  selection  of  particle  size  and 
concentration,  solution  properties,  electric  field  intensity  and  its  time 
protocol. 
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fig.  1:  Ajstefe'of  a  porous  substrate  of  thickness  c. 

fig.  2:  The  definition  of  the  coordinate  system  and  the  geometrical  model  for  particle 
penetration . 

fig.  3:  The  probability  density  Pd  vs.  the  downstream  nondimensional  distance  z/b. 
fig.  4:  The  probability  density  Pd  vs.  the  electrophoretic  Peclet  number, 
fig.  5:  The  probability  density  Pd  vs.  the  local  daa&ty  number, 
figs.  6-9:  The  cummulative  probability  Pc  vs.  the  downstream  nondimensional  distance 
z/b. 

fig.  10:  Penetration  depth  defined  as  the  downstream  value  distance  z/b  for  which  there  is 
90  percent  probability  that  the  particle  has  deposited  after  a  very  long  time  vs  the 
electrophoretic  Peclet  number  (the  local  deposition  number  used  as  a  parameter). 
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Abstract 


Zr02  coatings  were  deposited  on  graphite  and  titanium  from 
an  aqueous  solution  in  a  two-step  reaction  :  hydroxyl  ions 
(OH' )  were  generated  at  the  cathodic  substrate  by  reduction 
of  NO3 '  and  dissolved  O2  ,  and  then  reacted  with  zirconyl 
ions  in  the  solution  to  form  the  hydroxide  Zr(OH)-»  ,  which 
in  turn  decomposed  on  drying  to  yield  zirconia  (Zr02 )  - 
A  faradaic  efficiency  of  40-80%  was  found,  attributable  to 
reduction  reactions  that  do  not  produce  hydroxyl  ions  (such 
as  reduction  of  H‘  ),  as  well  as  to  formation  of  the 
hydroxide  at  sites  removed  from  the  cathodic  substrate.  The 
effects  of  current  density,  time  and  hydrodynamic  conditions 
on  coating  weight,  cell  voltage,  temperature  and  pH  of  the 
solution  were  studied. 

Introduction 

Ceramic  coatings  are  currently  of  much  interest  for 
applications  in  high-temperature  and  severe' y  corrosive 
environments.  A  variety  of  methods  are  available  for  their 
production,  including  plasma  spraying  f 1  - 2 1 ,  chemical  vapor 
deposition  and  sputtering  l1!,  and  sol -gel  processes  f3  si. 

Formation  of  ceramic  coatings  by  electrochemical  means  is 
a  relatively  new  technique  and  has  been  described  by  Switzer 
‘ 6 1  in  application  to  synthesis  of  ceramic  films  and 
powders.  It  presents  several  advantages  over  alternative 
coating  techniques  :  the  thickness  and  morphology  of  the 
deposit  can  be  controlled  by  the  electrochemical  parameters; 
relatively  uniform  deposits  are  obtainable,  on  complex 
shapes;  the  deposition  rate  is  higher  than  for  most  other 
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methods,  and  the  equipment  is  of  low  cost  l7l. 

Electrodeposition  has  been  used  extensively  for  many  years 
for  production  of  metallic  coatings.  Later,  oxides  and 
hydroxides  T6-91  were  deposited  by  a  redox  change  reaction 
induced  on  the  anode.  Formation  of  a  ceramic  oxide  through 
cathodic  reactions  was  demonstrated  only  recently  by  Switzer 
t  6 1  with  CeC>2  coatings. 

In  the  present  work,  deposition  of  zirconia  (ZrC>2  )  from  an 
aqueous  solution  of  zirconyl  nitrate  (Zr0(N03)2)  was 
studied.  The  sequence  of  reactions  leading  to  oxide 
formation  is  expected  to  be  the  following: 

a)  dissociation  of  the  zirconyl  salt  : 

Zr0(N03)2  — »  ZrO* 1  +  2  X03 

b)  hydrolysis  of  the  zirconyl  ion  : 

ZrO' '  +  H20  — >  Z r ( OH ) 2 ' 1 

c)  interaction  of  the  hydrated  cation  with  OH  ions 
generated  at  the  cathode  by  reduction  reactions 
described  further  on: 

Z r ( OH ) 2 1  1  +  2  OH-  — >  Zr(OH)i 

d)  dehydration  of  the  hydroxide  : 

Zr(OH)<  Zr02  +  2  H20 

The  cathodic  reactions  that  generate  OH  and  their  standard 
and  equilibrium  potentials  calculated  with  the  aid  of  the 
Nernst  equation  are  : 

e)  02  +  2  H20  +  4  e  — >  4  OH  E„  =  0.401  v 

E„q . =  1.0195  v 
(at  25  °C) 

02  +  4  H*  +  4  e  — >  2  H2Q  E„  =  1.299  v 

E„q  .  =  1.05  v 

£)  reduction  of  the  NO3  ions  : 

NO3-  +  H20  +  2  e  — »  NO2  +2  OH-  Eo  =0.01  v 

< — 

E,. q  .=  0.997  v 
(for  0.1  M  nitrate) 

g)  reduction  of  H2O  : 

2  H2O  +  2  e  H2  +  2  OH  Eo  =  0.828  v 

Eoq.— “0.165  v 

The  anodic  reaction  which  occurs  simultaneously  is  : 

H20  — »  1/2  O2  +  2  H'  +  2  e 

< - 

Thus,  the  hydroxide  is  formed  in  two  successive  steps  : 
electrochemical  generation  of  OH  ,  and  a  chemical  reaction 
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between  the  latter  and  the  zirconyl  cations.  The  deposit 
forms  on  the  cathode,  and  a  pla'tinum  electrode  serves  as 
anode . 

The  present  paper  describes  the  effects  of  the 
electrochemical  deposition  parameters  on  the  formation  rate 
of  the  Zr02  coatings. 

Experimental  Procedure 

The  deposits  were  obtained  on  two  different  substrates, 
prepared  in  the  form  of  2  x  10  x  20  mm  specimens  :  a)  porous 
graphite  (UC  grade  45)  and  b)  commercially  pure  titanium 
(ASTM  B-265-58T).  The  graphite  specimens  were  polished  with 
a  1000  grit  Sic  abrasive  paper,  and  the  Ti  plates  machined 
and  then  polished  with  180  grit  paper.  All  specimens  were 
then  rinsed  with  ethanol  in  an  ultrasonic  bath,  then  washed 
with  distilled  water  and  dried  in  air. 

The  electrolyte  used  was  a  0.1  M  aqueous  solution  of 
ZrO(NC>3  )2  •  nH2  0  ,  with  an  initial  pH  2.3  . 

A  power  supply,  Hipotronics  801-iA,  was  used,  and  cell 
voltage  and  current  were  measured  with  AVOmeters. 
Polarization  curves  were  obtained  with  a  PAR  model  270 
potentiostat .  Deposits  were  obtained  at  current  densities 
ranging  from  15  to  100  mA/cm2  ,  ana  durations  of  10-60 
minutes,  in  the  course  of  which  the  cell  voltage,  and  the  pH 
and  temperature  of  the  solution  were  measured.  Most 
experiments  were  performed  without  stirring,  some  with  it. 

The  deposits  were  dried  in  air  at  room  temperature. 
Coating  weights  were  determined  by  weighing  the  specimens 
before  and  after  deposition,  within  an  accuracy  of  +0.05  mg. 
Some  of  the  deposits  were  submitted  to  high-temperature 
treatments.  Those  on  graphite  were  treated  either  at  900 °C 
in  argon  or  at  400  and  600 °C  in  air,  while  those  on  Ti 
were  treated  at  780°C  in  air.  The  morphology  and  phase 
evolution  of  the  deposits  before  and  after  the  heat 
treatment  were  studied  with  the  aid  of  scanning  electron 
microscopy  and  x-ray  diffraction,  and  reported  in  H°l. 

Resul ts 

The  results  described  below  refer  only  to  the  graphite 
specimens.  Cathodic  polarization  curves  are  shown  in  Fig.l  . 
For  the  unstirred  solution,  two  polarization  regions  were 
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observed,  representing  reduction  reaction  associated  with 
activation-  and  concentration-type  polarization  regimes 
respectively.  Due  to  the  range  of  the  potentiostat  (limited 
to  -4  volt),  the  maximal  current  density  (c.d.)  achieved  was 
limited  to  55  mA/cm2  ,  which  prevented  detection  of  a  third 
cathodic  reaction  that  must  occur  during  deposition  at 
higher  c.d. 's. 

The  curve  obtained  for  the  stirred  solution  does  not  show 
domains  of  limited  c.d.,  but  the  regions  with  different 
slopes  indicate  three  cathodic  reactions,  which. will  be 
characterized  and  discussed  later  in  detail. 

The  dependence  of  coating  weight  on  the  c.d.  at  which  the 
coating  was  formed,  for  a  constant  deposition  time  of  60 
min,  is  shown  in  fig. 2  .  As  could  be  expected,  weight 
increases  with  current  density,  and  the  observed  three  modes 
of  dependence  apparently  reflect  the  three  cathodic 
reactions . 

The  dependence  of  coating  weight  on  deposition  time  for 
two  different  c.d.'s  is  shown  in  fig. 3  .  Coating  weight 
increases  with  time  in  a  decelerating  manner.  Stirring  is 
seen  to  cause  a  significant  reduction  in  the  deposition  rate 
(fig. 4),  by  facilitating  migration  of  OH  from  the  cathode 
into  the  bulk  of  the  solution,  so  that  the  hydroxide  formed 
does  not  accumulate  on  the  cathode.  The  effect  of  stirring 
is  weaker  at  high  c.d.'s  . 

At  high  c.d.'s  (above  75  mA/cm2),  "strings"  of  gas  bubbles 
were  observed  after  a  certain  time,  emerging  from  the 
surface  of  the  coated  substrate. 

The  deposition  yield  in  terms  of  faradaic  efficiency  is 
relatively  low  and  is  presented  in  fig. 5  as  the  actual 
coating  weight  compared  with  the  theoretical  weight, 
assuming  that  the  electric  charge  served  exclusively  for 
generation  of  OH-  ions. 

The  dependence  of  cell  voltage  on  deposition  time  for 
different  c.d.'s  is  shewn  in  fig. 6  .  The  voltage  is  seen  to 
increase  with  c.d.  for  a  constant  deposition  time.  At  higher 
current  densities,  the  voltage  stabilizes  at  shorter  Limes. 
Thi  behavior  is  related  to  the  mechanism  of  the  coating 
for.  ation  and  is  discussed  later.  .  1 

"  The  dependence  of  solution  temperature  on  deposit  and  time 
is  shown  in  fig. 7  ,  again  for  different  c.d.'s  .  It  is  seen 
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that  initially  (for  the  first  10  min  or  so),  the  temperature 
level  varies  within  a  fairly  narrow  range,  irrespective  of 
the  c.d.  ;  but  in  the  longer  run  it  rises  drastically  -  the 
higher  the  c.d.  ,  the  steeper  the  rise.  This  thermal  effect 
is  due  to  the  increase  in  electric  resistance,  following 
formation  of  the  deposit. 

The  time  pattern  of  the  pH  level  of  the  solution  is  shown 
in  fig. 8  .  Below  40  mA/cm2  the  pH  hardly  varied,  but  at 
higher  c.d.'s  significant  changes  were  observed  already  from 
the  outset.  The  explanation  may  be  that  the  pH  as  measured 
was  governed  by  the  H'  ions  generated  at  the  anode,  since 
the  OH*  ions  formed  at  the  cathode  were  used  up  in  deposit 
formation  without  contributing  to  the  resultant  pH  of  the 
bulk  solution  (see  also  below). 

Discussion 

In  keeping  with  the  electrochemical  potentials  of  the 
various  reduction  reactions,  the  first  reaction  in  the 
polarization  curve  is  reduction  of  O2  ,  with  a  relatively 
low  limiting  c.d.  because  of  the  low  concentration  of 
dissolved  O2  ;  the  second  reaction  is  reduction  of  the  NO3' 
ions,  and  the  third  (unreflected  in  the  polarization  curve) 
is  expected  to  be  reduction  of  H4  .  Thus,  while  the  first 
t.wo  reactions  produce  OH  ions,  the  third  results  in 
formation  of  H2  . 

With  respect  to  coating  weight  as  function  of  c.d. 
(fig. 2),  the  first  part  of  the  curve  (low  c.d.'s  of  0-15 
mA/cm2  )  is  characterized  by  a  relatively  moderate  (low 
deposition  rate)  slope  (0.7  mg . cm2  /mA),  in  turn  associated 
with  a  low  generation  rate  of  OH'  from  reduction  of  O2  (see 
fig.l).  The  steeper  slope  (of  the  order  of  2.0  mg.cm2/mA) 
in  the  second  part  (  15-75  mA/cm2 )  represents  a  higher  rate 
of  deposit  formation,  associated  with  a  high  rate  of  OH 
generation  from  the  reduction  of  NO3'  .  The  third  part  (  75- 
100  mA/cm2)  is  again  associated  with  a  relatively  moderate 
slope  (0.7  mg.cm2/mA),  due  in  this  case  to  onset  of 
reduction  of  H'  ions  to  form  H2 ,  a  reaction  that  does  not 
produce  OH"  ions,  hence  does  not  contribute  to  deposit 
formation.  As  c.d.  increases,  more  competitive  electron¬ 
consuming  reactions  are  activated.  , 

The  deviation  of  the  actual  coating  efficiency  from  the 
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theoretical  values  may  be  due  to  the  following  phenomena  : 

a)  Charge  transfer  reactions  that  do  not  produce  OH"  ions  . 

b)  Consumption  of  OH"  by  H*  ions  present  in  the  original 
solution  and  generated  at  the  anode. 

c)  Undeposited  Zr(OH)<  formed  in  the  bulk  of  the  solution. 

The  increase  in  cell  voltage  with  deposition  time  may  be 

explained  as  follows  : 

In  the  initial  stage  of  the  coating  process,  the  cathodic 
substrate  is  conductive  and  OH"  ions  are  generated  over  its 
entire  surface.  Due  to  the  porosity  of  the  substrate,  a 
porous  coating  is  initially  formed.  As  the  coating  process 
progresses,  the  deposit  becomes  denser  and  thicker,  and  its 
electrical  resistance  increases,  with  the  attendant  increase 
in  both  cell  voltage  and  solution  temperature.  The  change  in 
electric  resistance  is  responsible  in  part  for  the 
accentuated  limiting  c.d.'s,  and  the  relatively  high 
potentials  observed  in  the  cathodic  polarization  curve  of 
the  unstirred  solution  (fig..l).  After  a  certain  time,  the 
cell  voltage  reaches  maximum,  but  the  deposit  continues  to 
form  while  undergoing  local  breakdowns  manifested  in 
"strings"  of  gas  bubbles.  These  breakdowns  create  sites  of 
stronger  reduction  currents,  which  in  turn  enable  the 
deposition  process  to  continue,  without  further  increase  in 
the  cell  voltage.  The  deposit  thus  consists  of  two  layers 
(see  fig.l  in  [10]);  a  thin  dense  layer  formed  in  the  first 
stage  and  a  relatively  thicker  one  formed  in  the  second.  The 
morphology  of  the  deposit  is  in  agreement  with  the  proposed 
mechanism  for  formation  of  thick  coatings. 

The  temperature  increase  may  affect  the  diffusion  rate  of 
OH  away  from  the  cathode  and  that  of  Zr(OH)2''  toward  it. 
Thus,  it  would  be  desirable  to  control  the  temperature 
during  the  coating  experiments  by  stirring  the  solution; 
This,  however  would  interfere  with  deposit  accumulation  on 
the  cathode.  Coating  weight  of  nickel  hydroxide  was  reported 
to  be  unaffected  by  moderate  agitation  of  the 
concentrated  solution.  In  this  regard,  coating  weight  is 
expected  to  be  influenced  when  low-concentration  solutions 
and  severe  stirring  conditions  are  used.  It  should  be 
mentioned  that  the  pH  was  measured  in  the  bulk  of  the 
unstirred  solution;  local  readings  at  the  cathode  may  show 
an  increase  in  pH  due  to  the  OH  ions. 
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The  present  work  demonstrated  the  outstanding  potentiali¬ 
ties  of  the  electrolytic  method  for  producing  ceramic 
coatings,  with  relatively  easy  control  of  their  thickness. 
However,  it  is  necessary  to  control  the  drying  stage  in 
order  to  achieve  sound  and  crack-free  coatings.  Optimal 
conditions  for  the  formation  of  Zr02  coatings  are  low 
current  densities  and  short  coating  times,  in  view  of  the 
non-conductive  nature  of  the  oxide. 
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Fiaure  Captions 


Fig.  1  :  Cathodic  polarization  curves  for  graphite  in  0.1  M 
Zr0(N03 )2* nHiO  ,  at  2  mV/sec  scan,  stirred  vs. 
unstirred  solution. 

Fig.  2  :  Coating  weight  as  function  of  current  density  for 
constant  deposition  time  of  60  min. 

Fig.  3  :  Coating  weight  as  function  of  deposition  time  for 
different  current  densities. 

Fig.  4  :  Coating  weight  as  function  of  deposition  time, 
stirred  vs.  unstirred  solution. 

Fig.  5  :  Coating  weight  as  function  of  electric  charge. 

Fig.  6  :  Cell  voltage  as  function  of  deposition  time  for 
different  current  densities. 

Fig.  7  :  Solution  temperature  as  function  of  deposition 
time  for  different  current  densities. 

Fig.  8  :  pH  of  solution  as  function  of  deposition  time  for 
different  current  densities. 
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Fig.  1  :  Cathodic  polarization  curves  for  graphite  in  0.1  M 
Zr0(J?03  )2 •  J1H2G  ,  at  2  mV/sec  scan,  stirred  vs. 
unstirred  solution. 
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3  :  Coating  weight  as  function  of  deposition  time  for 
different  current  densities. 
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Fig.  4  :  Coating  weight  as  function  of  deposition  time, 
stirred  vs.  unstirred  solution. 
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Fig.  7  :  Solution  temperature  as  function  of  deposition 
time  for  different  current  densities. 
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Abstract 

Zirconia  coatings  were  formed  on  graphite  and  titanium 
substrates  by  electrodeposition.  The  as-deposited  coatings 
were  amorphous,  and  cracked  during  the  drying  stage. 
Sintering  of  the  coatings  on  the  graphite  substrates  caused 
crystallization,  with  formation  of  a  variety  of  mixtures  of 
the  tetragonal  and  monoclinic  Zr02  polymorphs,  while 
oxidation  of  the  Ti  substrate  and  reaction  with  zirconia 
during  the  heat-treatment  resulted  in  a  microstruture 
resembling  ceramic  composite  coatings. 

Introduction 

Ceramic  materials  are  finding  ever  wider  application  as 
thermal -barrier  ti-21,  wear-resistant  f  3  *  *  1  ,  and  corrosion- 
resistant  t 5 ~ 6  J  coatings. 

The  main  ceramic  coating  methods  fall  into  three 
categories  :  physical  deposition,  chemical  deposition  and 
reaction/conversion  f7l.  An  electrochemical  method  also 
employed  is  based  in  on  electrophoretic  process  in  which  the 
ceramic  powder  is  suspended  in  a  fluid,  and  deposited  under 
the  action  of  a  high-strength  electric  field  1 0 - x x 1  . 

Recently,  several  workers  used  an  electrochemical  method 
to  form  coatings  of  transition-metal  oxides  and  hydroxides 
fl2l,  including  thallium  oxide  and  cerium  oxide 
1155  on  metallic  and  semiconducting  substrates  directly  from 
an  aqueous  solution.  In  such  process,  the  microstructure  and 
composition  of  the  coating  can  be  controlled  by  varying  the 


electrochemical  parameters  and  the  solution  composition. 

The  present  paper  deals  with  the  microstructural  evolution 
of  the  electrolytically  formed  zirconia  coatings  on  graphite 
and  titanium  substrates,  as  described  in  part  I  of  this 
series  r 1 6 1 . 

Experimental  procedure 

The  substrates  chosen  were  2  x  10  x  20  mm  plates  of 
commercial  graphite  (48  %  porosity,  average  pore  size  of  60 
micron)  and  commercially  pure  titanium  (ASTM  #  3  265-58T) 

materials  which  differ  in  their  chemical  affinity  to 
zirconia.  The  graphite  was  polished  with  1000  grit  Sic 
abrasive  paper,  while  the  titanium  plates  were  surface- 
finished  by  machining.  The  plates  were  finally  rinsed  with 
ethanol  and  distilled  water  and  then  dried  in  air. 

The  electrolyte  was  a  0.1  M  aqueous  solution  of  zirconyl 
nitrate  ZrO(NC>3  )2‘ 11H2O,  with  an  initial  pH  adjusted  to  2.3  . 
A  total  surface  area  of  “4  [cm2]  of  the  substrate  was 
immersed  in  the  solution.  Deposition  was  effected  at  current 
densities  within  the  range  of  15-100  [mA/cm2]  and  10-60 
minute  durations.  For  a  more  detailed  description  of  the 
deposition  process  see  rn)  . 

The  deposits  formed  on  the  substrates  were  dried  in  air  at 
room  temperature.  Firing  was  carried  out  at  400,  600°C  in 
air,  and  900°C  in  argon  ,  for  1  hr  for  the  graphite 
substrates,  and  at  780 °C  for  1  hr  in  air  for  the  titanium. 

The  microstructure  of  the  zirconia  coatings  was 
characterized  during  the  different  stages  of  the  process, 
using  optical  and  scanning  electron  microscopy  (SEM)  (Model 
JSM-840).  The  phase  content  was  determined  by  X-ray 
diffraction  with  the  diffractometer  (Model  PW-1730)  operated 
at  40  KV  and  40  mA,  using  monochromatized  Cu  k-alpha 
radiation,  at  a  scanning  speed  of  0.4  [degree/min] . 

Resul ts 

The  as-deposited  ceramic  coatings  appeared  as  a 
transparent  viscous  fluid,  especially  at  high  current 
densities  and  long  coating  durations,  where  thick  deposits 
were  formed.  In-situ  microscopical  observations  of  the 
drying  process  showed  simultaneous  development  of  cracks 
within  the  coating  layer  on  graphite.  The  wet  coating 


deposit  was  characterized  by  a  high  content  of  bubbles. 
During  drying  some  of  the  bubbles  coalesced  to  form  larger 
ones,  while  others  were  expelled  from  the  deposit  and 
migrated  to  its  surface.  The  drying  process  was  thus 
accompanied  by  formation  of  discontinuities  within  the 
coating,  which  converted  into  deep  cracks  towards  the  end  of 
the  drying  process.  The,  resultant  coating  layer  had  a 
”cracked-mud"  appearance  (as  shown  in  figs.l  and  2  ) 
irrespective  of  the  substrate  material.  The  coating 
thickness  was  found  to  increase  with  both  current  density 
and  coating  duration.  Thicker  coatings  showed  larger  islands 
in  the  cracked-mud  morphology. 

Thick  coatings  (over  ~  20  micron)  on  the  graphite 
substrates  often  consisted  of  two  layers,  the  upper  layer 
being  thicker  and  composed  of  coarser  islands  (fig. lb). 
Circular  traces  of  the  drying  pools  were  also  observed, 
especially  at  the  centers  of  these  coarse  islands  (fig.lc), 
but  were  almost  absent  in  the  coatings  on  the  Ti  substrates 
(fig. 2b) . 

X-ray  diffraction  spectra  of  the  dry  as-deposited  coatings 
and  of  the  fired  coatings  are  shown  in  figs. 3  and  4  for  the 
graphite  and  Ti  substrates,  respectively.  The  coating  layers 
before  firing  have  the  typical  spectrum  of  amorphous 
zirconia,  in  which  a  broadened  amorphous-type  peak  appears 
around  the  angle  28  ~  30  degrees  (figs.  3a  and  4a)  irrespec¬ 
tive  of  the  substrate  material.  This  finding  is  consistent 
with  the  relatively  smooth  and  featureless  morphology  of  the 
as-deposited  coating  surfaces  (figs.l  and  2). 

Firing  induced  a  different  phase  evolution  and  morphology 
in  these  substrate/coating  systems,  and  the  resulting  fired 
coatings  were  opaque.  The  microstructure  of  the  coating  on 
the  graphite  substrate  was  characterized  by  fine  equiaxed 
submicron  particles  (fig. 5).  These  particles  were  identified 
as  a  mixture  of  tetragonal  and  monoclinic  polymorphs  of 
zirconia  (fig. 3b),  based  on  the  {111}  and  {400}-type 
reflections.  The  broadening  of  the  {lll}-type  reflections 
was  used  to  calculate  the  crystallite  size  of  the  polymorphs 
with  the  aid  of  the  Scherrer  equation  ti7]#  while  the  volume 
fractions  of  the  two  phases  were  determined  through  the 
'polymorph'  equation  as  corrected  by  Porter  and  Heuer  t1®!. 
The  crystallite  size  was  found  to  be  2  and  12  [nm]  for  the 


400 °C  and  600 °C  treatments  respectively,  with  the  polymorphs 

O 

undistinguishable.  By  contrast,  treatment  at  900  C  for  1  hr 
in  argon  resulted  in  distinct  crystallite  sizes  of  25  and  20 
[nm]  for  the  tetragonal  and-monoclinic  phases,  with  volume 
fractions  of  15%  and  85%  ,  respectively.  These  size  data 
indicate  that  the  SEM  image  (fig.  5)  represents  zirconia 
aggregates.  , 

The  zirconia  coating  on  the  Ti  substrates  exhibited  a 
totally  different  inner  microstructure.  Oxidation  of  the 
metallic  titanium  during  firing  resulted  in  the  original 
spaces  between  the  cracks  in  the  coating  layer  being  filled 
by  the  growing  titanium  oxide  (fig. 6).  Microchemical  EDS 
analyses  of  the  various  features  of  the  coating  (see  Table 
l),and  elemental  mapping  for  Ti  (fig.6d)  confirmed  the 
Ti-rich  nature  of  the  polycrystalline  phase  between  the 
zirconia  islands.  An  appropriate  X-ray  diffraction  spectrum 
from  this  specimen  (fig.  4b)  contained  reflections  of  Ti02 
in  the  form  of  rutile,  Zr02  with  the  tetragonal  and 
monoclinic  phases,  and  additional  reflections  coincident 
with  those  of  the  Zrc . 72Tio . 2aOo . 28  ,  Zro . s Tio . 5 Oo . 33 , 
Zro . 5 Tio . 5 Oo . 1 0  hexagonal  phases  f191.  The  volume  fractions 
of  the  tetragonal  and  monoclinic  phases  in  these  coatings 
were  calculated  to  be  37%  and  63%  ,  and  their  crystallite 
sizes  26  and  25  [nm] ,  respectively. 

The  surface  microstructure  was  more  complex  and  dependent 
on  the  coating  layer  thickness  prior  to  firing.  At  locations 
with  greater  coating  thicknesses  (i.e.  at  the  bottom  of  the 
machining  grooves),  titanium  oxide  growth  was  confirmed  to 
the  gaps  between  the  zirconia  islands  (fig. 6c).  However,  at 
locations  with  smaller  coating  thicknesses  (i.e.  top  of  the 
machining  grooves),  the  oxide  grew  to  extent  that  it 
partially  buried  the  zirconia  islands  (fig.6e).  In  both 
cases,  the  centers  of  the  islands  appeared  to  have  reacted 
with  the  titanium  substrate  (see  also  fig. 9),  and  were  most 
probably  converted  into  the  titanium  zirconate  phases.  EDS 
results  from  such  regions  (Table  I)  have  revealed 
compositions  similar  to  those  of  titanium  zirconate  phases, 
which  were  identified  by  XRD  (fig. 4b). 

In  many  cases,  traces  were  observed  on  the  titanium  oxide 
layer  parallel  to  the  periphery  of  the  zirconia  islands 
(arrowed  in  fig. 7).  These  traces  could  be  identified  as  the 


original  contours  of  the  zirconia  prior  to  firing.  By  their 
means,  the  firing  shrinkage  of  the  coating  was  estimated  to 
be  of  the  order  of  "20  %,  assuming  a  similar  percent 
shrinkage  for  the  third  dimension  (the  coating  thickness). 

The  diamond  Vickers  microhardness  of  the  coating  on  the  Ti 
substrates  has  an  average  value  of  520  [Kg/mm2 ]  .  Scratches 
were  applied  manually  to  the  coating  surface  with  a  diamond 
indentor,  in  order  to  characterize  qualitatively  the 
adherence  and  brittleness  of  the  composite  coating  layer.  A 
typical  trace  of  one  of  these  scratches  is  shown  in  fig. 8  ; 
it  exhibits  a  high  density  of  microcracks  perpendicular  to 
the  scratch  direction.  The  presence  of  these  microcracks  in 
the  smeared  region  indicates  heavy  deformation  of  the 
material  by  the  indentor.  Moreover,  some  regions  of  the 
coating  near  the  scratch  trace  were  debonded  from  the 
substrate  (fig. 9).  This  debonding  occurred  by  crack  propaga¬ 
tion  either  through  the  coating  or  at  the  substrate/coating 
interface.  In  the  latter  case,  the  fractured  surface  was 
composed  of  very  fine  dimples,  indicating  some  ductility  at 
the  interface. 

Discussion 

The  microstructural  and  morphological  evolution  within  the 
coating  layer  may  be  explained  in  relation  to  the  affinity 
of  the  substrate  material  to  chemical  reaction  with  the 
coating. 

The  electrochemical  deposition,  whereby  zirconium 
hydroxide  gel  is  precipitated  from  the  solution  at  ambient 
temperature  on  the  cathode,  leads  to  the  amorphous  character 
of  the  coating. 

The  bubbles  resulting  from  the  drying  process  were 
observed  only  in  coating  layers  formed  on  the  graphite 
substrates.  Consequently,  the  bubbles  within  the  wet 
coatings  may  be  related  to  the  entrapped  gases  within  that 
substrate,  due  to  its  porous  nature.  The  non-porosity  of  the 
Ti  substrates  prevented  such  bubble  formation.  Still,  the 
coatings  in  both  systems  showed  surface  microcracks  after 
the  drying  stage.  This  type  of  morphological  development  is 
attributable  to  non-uniform  contraction  of  the  wet  coating. 
Densification  of  the  latter  during  the  drying  process 
results  in  tensile  stresses  within  it,  the  level  of  these 


stresses  depending  upon  the  mismatch  in  thermal  expansion 
coefficient  and  elastic  constants  between  the  coating  and 
substrate  materials.  While  thicker  coatings  may  withstand 
biaxial  stresses  in  their  plane  better  than  thinner  ones, 
stresses  perpendicular  to  the  coating/substrate  interface 
are  larger  for  thicker  coatings,  and  may  cause  spalling.  The 
resultant  morphology  is  reflected  in  fewer  cracks  in  the 
thicker  coating  layers,  hence  in  formation  of  larger 
zirconia  islands. 

Firing  of  the  zirconia  coatings  on  graphite  substrates  at 
the  various  temperatures,  led  to  crystallization  of  the 
amorphous  coating.  According  to  the  Ellingham  diagram  £20], 
graphite  oxidizes  under  these  firing  conditions  to  form  CO 
gas,  but  no  reaction  with  the  zirconia  is  expected  to  occur. 
Low-temperature  firing  of  amorphous  zirconia  was  reported  to 
produce  low- temperature  polymorphs  £21-22],  The  pH  of  the 
mother  solution  was  also  reported  to  affect  the  ratio  of 
tetragonal/monoclinic  phase  contents  after  calcination  at 
400  to  600 °C  £23-24)  ,  at  low  pH  (less  than  about  5)  a 
mixture  of  these  low  temperature  polymorphs  is  expected  to 
be  formed,  as  was  confirmed  for  the  present  coatings. 
Formation  of  undoped  zirconia  via  solution  by  different 
techniques  often  results  in  a  mixture  of  the  polymorphs, 
their  relative  volume  fractions  depending  upon  crystallite 
size  £25-26],  The  crystallite  size  of  the  tetragonal 
polymorph  present  in  the  coating  (12-25  nm)  is  also 
consistent  with  the  reported  particle-size  range  within 
which  this  phase  is  stable  £21-221, 

By  contrast,  firing  of  the  coatings  on  Ti  substrates  has 
resulted  in  a  more  complex  morphology,  due  both  to  oxidation 
of  Ti  and  to  its  reaction  with  the  coating.  In  this  case  the 
firing  temperature  was  780 °C,  at  which  Ti  exhibits  an 
oxidation  regime  of  a  parabolic  type  [27-20],  Under  this 
regime,  oxidation  is  said  to  occur  by  diffusion  of  Ti  atoms 
towards  the  surface  of  the  growing  oxide  layer.  Thus,  while 
the  surface  of  the  Ti  substrate  at  the  cracks  is  exposed  to 
air  and  oxidation,  the  regions  below  the  zirconia  islands 
are  subjected  to  relatively  low  oxygen  partial  pressures, 
assuming  solid  contact  at  the  zirconia/Ti  interfaces. 
Accordingly,  at  these  interfaces,  oxidation  of  Ti  occurs  by 
interdiffusion  of  Ti,  Zr  and  oxygen.  The  titanium  zirconate 


phase,  with  its  low  oxygen  stoichiometry,  is  evidence' of 
this  reaction. 

It  was  found  that  the  reaction  between  the  Ti  and  the 
coating  was  limited  to  the  centers  of  the  zirconia  islands, 
where  significant  volume  contraction  was  observed.  This 
effect  is  attributable  to  differential  shrinkage  of  the 
islands  during  firing.  Due  to  the  temperature  gradient 
across  the  zirconia  layer,  the  islands  are  expected  to  bend 
upward  at  the  edges  during  sintering,  thus  losing  contact 
with  the  Ti  substrate  along  their  periphery.  Furthermore, 
oxidation  of  Ti  at  780 °C  occurs  mainly  by  diffusion  of  Ti 
rather  than  that  of  oxygen  through  the  growing  oxide  layer, 
a  circumstance  which  favors  formation  of  the  low-oxygen 
titanium  zirconate  phase. 

Finally,  the  qualitative  scratch  test  of  the  coating  and 
its  debonding  behavior  indicate  relatively  good  adhesion  to 
the  Ti  substrate.  In  this  regard,  the  chemical  reactions 
between  substrate,  coating  and  atmosphere  determine  the 
microstructural  evolution  within  the  coating.  The 
characteristic  microstruct ure  which  resulted  from  these 
reactions  suggests  possibilities  for  formation  of  composite 
ceramic  coatings. 
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Table  I  :  EDS  composition  results  for  ceramic 
coating  on  Ti  substrates.* 


Microstructural 

feature 

Chemical  composition  (mole%) 

Remarks*  * 

Zr 

Ti 

Zirconia 

48.0 

52.0 

PI  in 

island  center 

59.0 

41.0 

fig. 6e 

(reacted 

65.0 

35.0 

regions) 

69.0 

31.0 

76.0 

24.0 

Zirconia 

73.0 

27.0 

P2  in 

island  periphery 

85.0 

15.0 

fig. 6e 

(non-reacted 

87.0 

13.0 

regions) 

i 

Oxidized  regions 

9.0 

91.0 

P3  in 

between  zirconia 

0.0 

100.0 

fig. 6e 

islands 

0.0 

100.0 

*  Fired  at  780°C  for  1  hr  in  air. 

**  Examples  for  the  analyzed  regions. 
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Figure  Captions 

Fig.  1  :  SEM  images  of  dried  zirconia  coatings  on  graphite 
resemble  a  "cracked-mud"  morphology  both  for  (a) 
thin  coating  (25  mA/cm2  for  15  min)  and  (b)  thick 
coating  (25  mA/cm2  for  75  min).  Thick  coatings 
(i.e.  100  mA/cm2  for  60  min)  are  often  composed  of 
two  layers,  with  the  upper  layer  containing  drying 
pools  (arrowed  in  (c)).  (d)  Cross-section  shows 
penetration  of  the  coating  into  pores  of  the 
graphite. 

Fig.  2  :  SEM  images  of  dried  zirconia  coatings  on  Ti  showing 
the  "cracked-mud"  morphology.  This  featureless 
surface  morphology  is  characteristic  of  the 
amorphous  nature  of  the  coating,  (b)  At  higher 
magnifications,  no  drying  pools  are  visible  in 
coatings  on  Ti . 

Fig.  3  :  X-ray  diffraction  spectra  of  zirconia  coatings  on 
graphite,  (a)  Amorphous-type  broadened  peak  from 
the  as-deposited  coating  and  the  sharp  peak  from 
the  graphite  (gr)  substrate,  (b)  Reflections  from 
the  tetragonal  (t)  and  monoclinic  (m)  polymorphs 
are  characteristic  of  the  Zr02  coatings,  sintered 
at  600 °C  for  1  hr. 

Fig.  4  :  X-ray  diffraction  spectra  of  zirconia  coatings  on 
Ti  .  (a)  Absence  of  sharp  reflections  from  a 
deposited  coating  indicates  on  its  amorphous 
nature,  (b)  Typical  reflections  from  alpha-Ti,  TiC>2 
(rutile-R),  Zr02  ( tetragonal -t  and  monoclinic-m  ) 
and  various  zirconium  titanates  (ZT)  are  seen  for 
coatings  sintered  at  780°C  for  1  hr. 

Fig.  5  :  SEM  image  of  sintered  (600°C/  1  hr)  zirconia 
coating  on  graphite.  The  submicron  particles 
represent  aggregates  constituting  a  mixture  of 
tetragonal  and  monoclinic  polymorphs  of  zirconia. 
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Fig.  6 


Fig.  7 


Fig.  8 


Fig.  9 


:  SEM  images  of  sintered  zirconia  coatings  on  Ti, 
showing  formation  of  Ti02  at  the  cracks  between  the 
zirconia  islands,  (a)  At  low  and  (b)  higher 
magnifications,  the  topography  of  the  machining 
gr.oves  is  visible,  (c)  At  the  bottom  of  grooves 
"BG"  ,  the  Ti02  crystals  grew  only  between  the 
zirconia  islands,  (d)  Ti  elemental  mapping  of  (c). 
(e)  At  the  top  of  grooves  "TG",  TiC>2  crystals  grow 
and  bury  the  zirconia  islands  (see  text  for  further 
details) . 

:  SEM  image  showing  the  line  traces  of  the  zirconia 
island  contours  prior  to  (arrowed)  and  after 
sintering,  from  which  the  sintering  shrinkage  may 
be  deduced. 

:  SEM  image  of  the  scratch  trace,  showing  high 
density  of  microcracks  at  the  smeared  regions, 
perpendicular  to  the  scratch  direction. 

:  SEM  images  showing:  (a)  Typical  coating  cross- 
section  where  fine  Ti02  crystals  grew  around  the 
zirconia  islands,  accompanied  by  shrinkage  and 
bending,  (b)  Debonding  at  the  coating/substrate 
interface  results  in  dimpled  fracture  surfaces. 


Fig.  1 


:  SEM  images  of  dried  zirconia  coatings  on  graphite 
resemble  a  "crackad-mud"  morphology  both  for  (a) 
thin  coating  (25  mA/cm2  for  15  min)  and  (b)  thick 
coating  (25  mA/cm2  for  75  min).  Thick  coatings 
(i.e.  100  mA/cm2  for  60  min)  are  often  composed  of 
two  layers,  and  the  upper  layer  containing  drying 
pools  (arrowed  in  (c)).  (d)  Cross-section  shows 
penetration  of  the  coating  into  pores  of  the 
graphite. 


£:ig.  2  :  S£M  images  of  dried  zirconia  coatings  on  Ti  showing 
the  "cracked-mud"  morphology.  This  featureless 
surface  morphology  is  characteristic  of  the 
amorphous  nature  of  the  coating,  (b)  At  higher 
magnifications,  no  drying  pools  are  visible  in 
coatings  on  Ti . 


6 


CM 


2  *  THETA  (degree ) 


X-ray  diffraction  spectra  of  zirconia  coatings  on 
graphite,  (a)  amorphous-type  broadened  peak  from 
the  as-deposited  coating. and  the  sharp  peak  from 
the  graphite  (gr)  substrate,  (b)  Reflections  from 
the  tetragonal  (t)  and  monoclinic  (m)  polymorphs 
are  characteristic  of  the  Zr02  coatings,  sintered 
at  600°C  for  1  hr. 


Fig.  5  :  SEM  image  of  sintered  (600 °C/  1  hr)  zirconia 
coating  on  graphite.  The  submicron  particles 
represent  aggregates  constituting  a  mixture  of 
tetragonal  and  monoclinic  polymorphs  of  zirconia. 


Fig.  6  :  SEM  images  of  sintered  zirconia  coatings  on  Ti, 

showing  formation  of  Ti02  at  the  cracks  between  the 
zirconia  islands,  (a)  At  low  and  (b)  higher 
magnifications,  the  topography  of  the  machining 
grooves  is  visible,  (c)  At  the  bottom  of  grooves 
"BG"  ,  the  Tio2  crystals  grew  only  between  the 
zirconia  islands,  (d)  Ti  elemental  mapping  of  (c). 
(e)  At  the  top  of  grooves  "TG",  Ti02  crystals  grow 

and  bury  the  zirconia  islands  (see  text  for  further 
details) . 


Fig.  7  :  SEM  image  showing  the  line  traces  of  the  zirconia 
island  contours  prior  to  (arrowed)  and  after 
sintering,  from  which  the  sintering  shrinkage  may 
be  deduced. 


Fig.  8  :  SEM  image  of  the  scratch  trace,  showing  high 
density  of  microcracks  at  the  smeared  regions, 
perpendicular  to  the  scratch  direction. 
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